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Isocurvature bounds on axions revisited 
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The axion is one of the best motivated candidates for particle dark matter. We study and update 
the constraints imposed by the recent CMB and LSS experiments on the mass of axions produced 
by the misalignment mechanism, as a function of both the inflationary scale and the reheating 
temperature. Under some particular although not unconventional assumptions, the axionic field 
induces too large isocurvature perturbations. Specifically, for inflation taking place at intermediate 
energy scales, we derive some restrictive limits which can only be evaded by assuming an efficient 
reheating mechanism, with T^h > 10^^ GeV. Chaotic infiation with a quadratic potential is still 
compatible with the axion scenario, provided that the Peccei-Quinn scale fa is close to 10^" or lO'^^ 
GeV. Isocurvature bounds eliminate the possibility of a larger fa and a small misalignment angle. 
We find that isocurvature constraints on the axion scenario must be taken into account whenever 
the scale of infiation is above 10^^ GeV; below this scale, axionic isocurvature modes are too small 
to be probed by current observations. 

PACS numbers: 



I. INTRODUCTION 

The fact that the strong sector of the Standard Model 
conserves the discrete symmetries P and CP while the 
electroweak sector doesn't, also known as the strong CP 
problem, is considered a serious puzzle for modern par- 
ticle physics The most elegant and compelling so- 
lution to this problem was proposed in 1977 by Peccci 
and Quinn 0] with the introduction of a new C/(1)pq 
global symmetry at high energies. The Peccei-Quinn 
(PQ) scalar field ip 
form 



As the universe expands, its energy decreases to about 
Aqcd ~ 200 MeV, and non-perturbative instanton ef- 
fects tilt the previously flat axion potential, explicitly 
breaking the residual symmetry Q: 



(3) 



p/v^e*® has a potential of the 



The axion field acquires a mass about the minimum of 
the potential that depends on the temperature in the 
vicinity of T Aqcd as 0, 13 



P 



7V2 J 



VpQ~lmlp' + ^p\ (1) 



(4) 



where N is the number of degenerate QCD vacua associ- 
ated with the color anomaly of the PQ symmetry. Spon- 
taneous symmetry breaking (SSB) occurs when the en- 
ergy density of the universe falls below VpQ and the field 
acquires a vacuum expectation value (VEV) p = fa/N. 
The axion 0] is then the Goldstone boson of the broken 
PQ symmetry. 
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where C is a model-dependent factor calculated in Refs. 
0, 9 to be of the order of C ~ 0.018. The mass 
finally reaches its asymptotic zero-temperature value 
rua = ma(0): 

ma = — —=6.2pcV[ . 5 

l + Z fa \ fa J 

Here z — rriu/md ~ 0.56 is the mass ratio of up to down 
quarks, while TOtt and /^r are respectively the pion mass 
and decay constant. The field equation of motion is 



Note, however, that SSB is effective only when the typi- 
cal fluctuations on 6p are smaller than fa/N. If either T 
or iJinf arc of order fa/N at reheating or during inflation, 
thermal [J or quantum fluctuations (respectively) will 
modify the effective potential and restore the PQ sym- 
metry, with a coherence length comparable to the Hubble 
scale at that time. 



0, 



(6) 



where V'{a) ~ dV/da, is the comoving laplacian, 
and R is the scale factor of the universe. If the axion 
field is initially displaced from the minimum of the po- 
tential when it acquires a mass, it starts oscillating in the 
potential described by Eqs. Q, Q, ©. 



The axion coupling to the rest of matter is inversely 
proportional to the symmetry breaking scale : 



contribute to the current dark matter component of the 
universe with a relic density [24| : 



Ant =57 :^E-B + tgfmf——f-f5f, (7) 

71" Ja Ja 

where a is the electromagnetic coupling constant, E and 
B are the electric and magnetic fields, / is a generic 
fermion and for each species gi is a model-dependent 
coefficient of order one. The main distinction between 
the various axion models comes from their coupling with 
electrons: for the "hadronic" models, such as the KSVZ 
model one has = 0, while the tree-level couphng 
does not vanish for the non-hadronic DFSZ models p^ . 
The other couplings are of the same order. For example, 
(7-y = —0.36 in the DFSZ model, while = 0.97 in the 
KSVZ model. 

In the currently accepted invisible axion model, the 
scale fa is in principle arbitrary, well above the elec- 
troweak scale so that the axion coupling to matter is 
weak enough to pass undetected, for the moment. There 
are at present several experiments searching for the axion 
in the laboratory, hke ADMX [H and CAST O, which 
have recently reported bounds on the axion coupling to 
matter [l^. 

Since in the large fa limit the axion remains effectively 
decoupled from other species, its fluctuations during in- 
flation could induce isocurvature perturbations in the 
CMB anisotropy spectrum which would in principle be 
observable today 0, [13, M, El, M, In this 

paper, we review the consequences of up-to-date bounds 
on the allowed isocurvature fraction and CDM density 
for the axion mass. 

The paper is organized as follows: in Section |TT] we 
review the different ways in which axions could be pro- 
duced, with a particular emphasis on the misalignment 
angle mechanism. In Section IIIII we study the induced 
isocurvature perturbations. In Section|TV]we present the 
various constraints bounding the axionic window. Fi- 
nally, our results are summarized and discussed in Sec- 
tion E] 



II. PRODUCTION MECHANISMS 

Axions arc produced in the early Universe by various 
mechanisms. Any combination of them could be the one 
responsible for the present axion abundance. We will 
briefly describe here the different production channels. 
For detailed reviews see Refs. [H, Q. 



A. Thermal production 

If the coupling of axions to other species is strong 
enough (i.e. fa low enough), axions may be produced 
thermally in the early universe and could significantly 



with = Pal Pc = ^T^Gpa/'iHQ and .g*,F is the number 
of degrees of freedom at the temperature at which the 
axions decouple from the plasma. The current WMAP 
bound on the dark matter density [2^ 

17edm/i' -0.1I2tO;™3, (9) 

together with Eq. ([5]), imposes a bound on the axion 
mass rua < 14.5 eV. As we will see in Section flVBl this 
bound is overseeded by astrophysical data which forbid 
a mass range of 0.01 eV < ma < 200 keV for the DFSZ 
axion P, Q , implying that the relic density of thermally 
produced DFSZ axions is completely negligible. 

Hadronic axions are not so tigthly constrained by as- 
trophysical data because they do not take part (at tree 
level) in the processes that cause the anomalous energy 
loses in stars such as 7 -I- e~ — > e~ 4- a or the Pri- 
makoff effect However, using the fact that thermally 
produced axions behave as warm dark matter, the au- 
thors of Ref. [2^ derived a model-independent bound 
nia < 1.05 eV showing that in any case i^^^"'^ ^ f^cdm- 
We will therefore ignore from now on the thermal axion 
contribution to dark matter, and work in the limit in 
which axions arc completely decoupled from the rest of 
matter. 



B. Production via cosmic strings 

We already mentioned that the PQ symmetry could be 
restored at high energy. After each symmetry restoration 
phase, a SSB can produce a population of axionic cosmic 
strings at the following epochs: 

• after inflation: if the scale fa/N is below the re- 
heating temperature of the universe, the PQ sym- 
metry is restored at reheating by thermal fluctu- 
ations. Axionic cosmic strings are produced later, 
when the temperature drops below fa/N These 
strings typically decay into axion particles before 
dominating the energy density of the universe. The 
axions produced this way are relativistic until the 
QCD transition, where they acquire a mass and be- 
come non-relativistic. Eventually these axions may 
come to dominate the energy density after equal- 
ity, in the form of cold dark matter. Estimates 
of their present energy density vary depending on 
the fraction of axions radiated by long strings ver- 
sus string loops. Three groups have studied this 
issue and found agreement within an order of mag- 
nitude [l,[23,[2i]. 



where Aqcd = 3 is a "fudge factor" which takes 
into account all the uncertainties in the QCD phase 
transition. Similar bounds were found in [29| fol- 
lowing a very different approach. 

• at the end of inflation: the PQ symmetry is re- 
stored during inflation whenever the typical ampli- 
tude of quantum fluctuations Hwif/'^T^ exceeds the 
symmetry breaking scale fa/N [5|. If the inequality 
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(11) 



holds throughout inflation, comic strings will be 
produced at the very end of this stage. The mech- 
anism by which cosmic strings are produced at the 
end of inflation is very different from that of a ther- 
mal phase transition and could affect the number of 
infinite strings, and thus the approach to the scal- 
ing limit and the relic density of axions. A detailed 
analysis of the axionic string production after in- 
flation and the corresponding estimate of the relic 
density of axions lies somewhat out of the scope 
of this paper and is left for future work. In the 
next sections and in Figs. [H [31 we will tentatively 
assume that the scaling limit is approached and 
that the corresponding relic density of axion is still 
given by Eq. PU)) . but the actual density could be 
significantly smaller. In any case, we will see that 
this bound is overseeded by the one derived from 
isocurvature modes. 

• during inflation: the PQ SSB could occur during 
inflation, when Hinf/in falls below fa/N. In this 
case, one expects axionic strings to be diluted dur- 
ing the remaining inflationary stage, and the relic 
density will be suppressed by an additional factor, 
exp(A^SSB), where A^ssb is the number of e-folds 
between PQ symmetry breaking and the end of in- 
flation. As a consequence, this density should be 
negligible today, unless A^ssb is fine-tuned to small 
values by assuming that fa/N is very close to the 
Hubble rate at the end of inflation. 

• before inflation: if, for instance, the PQ symmetry 
is restored at very high energy and breaks before 
inflation is turned on e.g. at low scales, the axionic 
strings produced in that way, as well as the possible 
axions into which they may have decayed, will be 
diluted by inflation and can be safely neglected. 

In addition, if iV > 1, axionic domain walls could be 
generated during the explicit symmetry breaking occur- 
ring around the QCD scale. Their energy density would 
end up as the dominant component of the universe [s^l 
unless one of the possible vacua turned out to be slightly 
non-degenerate. Since this would require a considerable 
amount of fine tunning, we will assume in the rest of this 
paper that = 1. A detailed discussion of this issue is 
given in Ref. [3l| . 



C. Generation via misalignment angle 

When the PQ symmetry is explicitly broken by instan- 
ton effects, the phase of the field 8 may not be at the 
minimum of its potential. As explained above, a (or Q) 
is a massless field during inflation and thus it fluctuates 
quantum-mechanically. If the typical amplitude of quan- 
tum fluctuations is large enough, Q could take different 
values in different points of our observable universe after 
inflation, with a flat probability distribution in the range 
[— TT, tt]; otherwise, it could remain nearly homogeneous. 
In both cases, at the time of the QCD transition, the (lo- 
cal or global) value of the misalignment angle 0qcd can 
differ from zero, leading to the sudden appearance of a 
potential energy m^/^(l — cosGqcd)- After the explicit 
symmetry breaking, the axion energy reads 
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for smaU 9.(12) 



Actually, the gradient energy can be safely neglected in 
Eq. (fT2|) . Indeed, even in the case in which the axion 
is maximally inhomogeneous, i.e. when the phase Q is 
equally distributed in the range [— tt, tt] in our observable 
universe at the end of inflation, it is straightforward to 
show that at any later time the coherence length (the 
physical size of the "homogeneity patches" for a) is al- 
ways of the order of the Hubble radius. This can be 
checked e.g. by solving the equation of motion ([6|) in 
Fourier space. As a consequence, and recalling that a is 
defined in the range 7r/a, 7r/a], the typical size of the 
gradient {S/a/R) is given by faH. So, the gradient energy 
scales as oc during radiation domination, and at 
the time of the QCD transition it is at most of the order 
of ifaHQco)^- A quick estimate gives Hqcb ~ 10~"eV, 
while in the following we will always consider values of 
the axion mass much larger than this. So, when the axion 
mass is "switched on" , the gradient energy is negligible 
with respect to the potential energy V ^ (fa'ma)'^- 

When ma{T) grows suddenly to values much bigger 
than iJqcD , the field quickly rolls down towards the mini- 
mum of the potential and starts oscillating with an en- 
ergy density given in first approximation by 
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Pa = {^f![e^ 



^^e^)) = i/>^(e|cD) (^) 

(13) 

However, for a precise estimate, it is necessary to take 
into account the time dependence of the mass. This was 
done in Ref. |3^], which found an extra factor fc ~ 1.44 in 
Eq. (Uni). Finally, using the conservation of the number 
density = Pa/maiT) in a comoving volume, one finds 
a relic axion density 



o h^^roA -5/12,^2 ^ /200McVy 
n^h ~ 7.245, / (0qcd) ^^^^ j 
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(14) 



where g^,i ~ 61.75 is the number of relativistic degrees 
of freedom in the universe at the temperature when the 
axion starts oscillating. Note that in Eq. (|14p. the spatial 
average (0qcd) °f ^he initial misalignment squared angle 
over the observable universe is not given by any field theo- 
retical reasoning, but by considerations on the stochastic 
behavior of the field during inflation. Using the Fokker- 
Planck equation, it is easy to show that 0, HH, 113 

mn-eQcu)Y^'-^x^AN^, (15) 

where Qi^ is the average value of the phase in our in- 
flationary patch at the time of Hubble crossing for the 
observable universe (A'obs e-folds before the end of infla- 
tion), and AiVqcD — 30 is the number of e-folds between 
TVobs and the time of Hubble exit for the comoving scale 
which re-enters the horizon when H = Hqcb ■ 

At this point, we see that two situations can occur. 
First, if fa 3> Hini, the right-hand-side in equation pS]) 
can be much smaller than one at the end of inflation; 
then, the axion field is essentially homogeneous, and the 
background value Oqcd — ©in in our universe is random 
but unique. Second, if fa < -ffinf , tbe right-hand-side can 
be of order one or larger, which means that the Brownian 
diffusion of the axion is complete, and the misalignment 
angle at the QCD scale is randomly distributed with a flat 
probability distribution in the range [— tt, tt]. Note that 
in this case, the quantum perturbations of the radial part 
of the PQ field are also large during inflation. In both 
cases, the mean energy density of the axion around Aqcd 
is proportional to 

(0QCd) = ^ £ c^da = y , (16) 

where the average should be understood as holding over 
many realizations of the universe in the case of a nearly 
homogeneous ©qcd, or over our present Hubble radius 
in the case of complete diffusion. 

Up until this point we have ignored the anharmonic 
corrections that could arise from the possibly large value 
of Oqcd- The calculations have been made using the 
approximation (1 — cos Q) ~ ^0^, which is obviously not 
valid for large angles. Including anharmonic corrections, 
one finds an enhancement factor 1.2 in Eq. ((TH)) . 



III. ISOCURVATURE PERTURBATIONS FROM 
AXION FLUCTUATIONS 

If the PQ symmetry is spontaneously broken during 
inflation, while the scale of inflation is much higher than 
that of the quark-hadron transition, the flat direction as- 
sociated with the massless Nambu-Goldstone boson will 
be sensitive to de Sitter quantum fluctuations. Indeed, 
quantum fluctuations are imprinted into every massless 
scalar field present during inflation, with a nearly scale 



invariant spectrum. 

If the scale of inflation is high enough, i7inf/27r > /a, 
it is possible that quantum fluctuations of the radial 
part of the PQ field restores the symmetry This 
symmetry restoration could have very different impli- 
cations for cosmological perturbations than a possible 
thermal symmetry restoration taking place after infla- 
tion. Indeed, the effective mass-squared V"{p) at the 
false vacuum p = is much smaller than the Hubble rate, 
V"[p) = f/487r2 < i/? f . So, the PQ field behaves like 
a light complex field during inflation. The symmetry is 
restored in the sense that an average over a scale much 
larger than the coherence length [t) of the field would 
give {tp) = 0. However, A^(t) is of the same order as the 
Hubble radius c/H{t) at a given time. In comoving space, 
the coherence length decreases by a huge factor dur- 
ing inflation, and the evolution of the field can be seen 
as a stochastic process of fragmentation into smaller and 
smaller homogeneity patches. But at the time when our 
observable universe crosses the Hubble scale, the PQ field 
is still nearly homogeneous inside our patch. Its quan- 
tum fluctuations become frozen beyond the horizon, and 
could thus leave a long wave perturbation which would 
still be described by Eq. (fT7|) . A detailed proof of this 
highly non-linear process requires a lattice simulation, 
whose analysis we leave for a future publication. In this 
paper, we will conservatively assume that isocurvature 
perturbations are erased when i?inf /2tt > fa. as has been 
assumed so far in the rest of the literature @ . 

The axion field perturbations 6a do not perturb the 
total energy density, first because the potential en- 
ergy is exactly zero, and second because, as explained 
above, the gradient energy of the axion cannot exceed 
^ (faHini)^; for fa ^ Mp this is much smaller than 
the total energy density (3/87r)Mp_ffj-^f . Since the to- 
tal energy density is unperturbed by these fluctuations 
during inflation, they are of isocurvature type, and man- 
ifest themselves as fluctuations in the number density of 
axions [ll,[i3,[ll,lil,[23l, 

5(^)^0. (18) 

In the absence of thermal symmetry restoration after in- 
flation, i.e. if the temperature of the plasma does not 
reach fa, the axion does not couple significantly to or- 
dinary matter. Its perturbations remain truly isocurva- 
ture [33| and may contribute as such to the temperature 
anisotropies of the CMB [sl, HI, . This isocurvature 
mode is expected to be completely uncorrelated to the 
usual adiabatic mode seeded by the quantum fluctuations 
of the inflaton. 

Let us assume that the Universe contains photons (7), 
approximately massless neutrinos (z^), baryons (b), ax- 
ions (a), ordinary CDM such as ncutralinos (x) and a 



cosmological constant. In the following, the subscript 
cdm will denote the total cold dark matter component, 
so that ficdm = ^a. + ^x- For the isocurvature mode, the 
perturbation evolution starts from the initial condition 
|(5^ = jSi, = 5h = 5x — and (Sa = S^, where 5a is the 
gauge invariant entropy perturbation 
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(19) 



(indeed, after the QCD transition, the axion fluid is non- 
relativistic with pa = ^a^-a, so (5a = (<5'^a)/'T-a; further- 
more, the fact that (5a 3> S-y implies Snu/ria 3> AST /T and 
5a = Srii^/nii = (5a)- It is equivalent to consider the per- 
turbations of a single cold dark matter fluid, obeying now 
to the initial condition (5cdm — -Ra'5a + (1 ~ Ra.)Sx ~ RaSa- 
If we compare with the initial condition for a usual 
"Cold Dark matter Isocurvature" (GDI) model, given by 
<5cdm = i^cdm, we sec that the axionic isocurvature solu- 
tion is equivalent to the GDI solution with 5cdm = RaSa- 
In other words, an axionic model with axionic fraction 
i?a = i^a/i^cdm, initial curvature spectrum (7?.^) and ini- 
tial entropy spectrum (5^) is strictly equivalent to a 
mixed adiabatic+GDI model with the same curvature 
spectrum and {S^^J = RUS^). 

Let us now relate the curvature and entropy power 
spectrum to the quantum fluctuations of the inflaton and 
axion field during inflation. For the adiabatic mode, it is 
well-known that the curvature power spectrum reads 



(20) 



where e is the first inflationary slow-roll parameter [37| 
and the subscript k indicates that quantities are evalu- 
ated during inflation, when k = aH. In first approxima- 
tion this spectrum is a power-law with a tilt riad depend- 
ing also on the second slow-roll parameter rj [s^l , 



Tiad = 1 - 6efe + 2rik 



(21) 



For the isocurvature mode, using the axion perturbation 
spectrum of Eq. PT|) , we obtain 
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(22) 

This power spectrum can be approximated by a power- 
law with a tilt 



1 - 2efc 



(23) 



which is related to the tilt rit of tensor perturbations, 
"iso = 1 + n*. 

The relative amplitude of isocurvature perturbations 
at a given pivot scale in adiabatic-|-GDI models is often 
parametrized as [H, [H, 110, Hlj 



(|5cdm(fc)|^ 



(|5edm(fc)P) + (|7^(fc)|2) 



(24) 



Since the axionic model is equivalent to an adia- 
batic-|-GDI model, we can still use the same parametriza- 
tion. The parameter a is related to fundamental param- 
eters by 
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(25) 



Rl{\SaikW) + {\n{kW, 'Wa\-QCD/ 

where in the last equality we assumed a <^ 1. 

IV. THE AXIONIC WINDOW 

We now describe the generic constraints on fa coming 
from different cosmological and astrophysical considera- 
tions. Some of them are generic and apply to every pro- 
duction mechanism or inflationary scenario while some 
others are rather model dependent. More precisely, we 
will enumerate the different bounds on the axion param- 
eter space (Minf, fa) associated with the misalignment 
mechanism of axion production during inflation. Mint be- 
ing the energy scale during inflation 



(26) 



and A/p = G^^/^ is the Planck mass. The explicit impli- 
cations of these bounds are presented on Figs. [U [31 



A. The scale of inflation 

The non-detection of tensor modes imposes a model- 
independent constraint on the inflationary scale. 



Minf < 3 X 10^'' GeV, 



(27) 



which follows from the current limit on the tensor to 
scalar ratio [25j 

/4 X lO'^iJinfA^ ^ , , 

r - — — < 0.3 at 95% c.l. . (28) 

V Mp / 

In Figs. [21 [21 this constraint corresponds to the upper 
hatched forbidden region. 



B. Supernova 1987A bounds 

The observed neutrino luminosity from supernova 
1987A imposes a bound on the axion luminosity that is 
saturated for 10"^ eV < iria < 2eV @. Other astrophys- 
ical and laboratory searches rule out an axion heavier 
than 1 eV (see [H, Q for a detailed discussion). There- 
fore, we have: 



nia < 10-^ - 10"'^ eV, 



(29) 



or equivalently, > 10^ - 10^° GeV. The forbidden re- 
gion is (blue-) shaded in Figs. [21 [31 



C. Axionic cosmic strings production 

As mentioned before, in the case of symmetry restora- 
tion at high energy, axions can be produced from the 
decay of cosmic strings, and we should impose a bound 
on their reUc density 

toT < COcdm , (30) 

where Wi stands for flih^ and ri|*'' is taken from Eq. (fT(I|l . 
Taking flcdmh'^ < 0.123 (see Eq.(|43l)) we get 

fa < 1.25 X 10"GeV . (31) 

This incquahty must be imposed in two cases: fa > 
(ffinf/27r), corresponding to symmetry restoration during 
inflation^ 0], and fa > Trh, corresponding to symmetry 
restoration after reheating [j| . ActuaUy, it is worth men- 
tioning that thermal corrections induce a positive mass- 
squared term m^g = so the precise condition for 
symmetry restoration is 

T^J12 » A/2 . (32) 

Therefore, the usual statement that symmetry is restored 
whenever fa > Tih assumes that the self-coupling con- 
stant A is of the order of 0.1. We will proceed with this 
assumption, but one should keep in mind that the exact 
condition is modcl-dcpcndcnt. 

The actual reheating temperature Tih is still unknown. 
This is the temperature at which the inflaton decays, 
once its half life has been exceeded by the age of the uni- 
verse, that is, when ~ F. Most of the thermal energy 
comes from perturbative decays of the inflaton and, as- 
suming that the decay products are strongly interacting 
at high energies, we can estimate the reheating temper- 
ature as 

^ 0.1 Vo7p ~ 0.02 /leff Vm Mp < 2 x lO" GeV , 

(33) 

where F = h'^gin/STr is the inflaton decay rate, which 
is typically proportional to the inflaton mass, m, with 
^eff < 10"'^, in order to prevent radiative corrections 
from spoiling the required flatness of the inflaton poten- 
tial [431. This estimate shows the generic inefficiency of 
reheating after inflation, where the scale of inflation could 
be of order 10^'^ GcV and the reheating temperature ends 
being many orders of magnitude lower. For instance, for 
Starobinsky type inflation, the weak gravitational cou- 
plings give a reheating temperature of order Ti± ~ 10^ 
GeV, while in chaotic inflation models, typical values are 



Note that for simplicity, wc impose this condition as if was 
constant at least during the observable e-folds of inflation (typi- 
cally, the last sixty e-folds). In principle, the amplitude of quan- 
tum fluctuations {Hi^{/2iT) could fall below fa precisely during 
the observable e-folds, see e.g. [T7| . but we will ignore this pos- 
sibility here. 



of order 10^" — 10^^ GcV. On the other hand, in certain 
low scale inflationary models, such as hybrid inflation, 
the efficiency of reheating can be significant because the 
rate of expansion at the end of inflation is much smaller 
than any other mass scale and the inflaton decays before 
the universe has time to expand, therefore all the inflaton 
energy density gets converted into radiation. 

We can parametrize the effect on the rate of expan- 
sion by introducing an efficiency parameter, ees, such 
that i/rh = ^cS Hcnd- Valucs of e range from 10"-'^^ for 
Starobinsky inflation, to order one for very low scale in- 
flation. If the reheating temperature is higher than fa 
it could eventually lead to a restoration of the PQ sym- 
metry. The subsequent spontaneous symmetry breaking 
would generate axionic cosmic strings that would not be 
diluted away by inflation. 

In summary, if the PQ symmetry is restored by thermal 
fluctuations after reheating, i.e. 

T,h = 0.l^H,i,Mp = O.lv/ecffi^infMp > fa , (34) 

then we must impose the condition ([30]) . In Figs. [21 [3l we 
have distinguished two cases: one in which the process of 
reheating the universe is very inefficient (toff < 10"^^), 
and there is no thermal restoration of the PQ symmetry 
after inflation, and another one in which eoff = 10~^ so 
that the symmetry might be restored. In both cases the 
constraint coming from string production and decay cor- 
responds to the triangular (red-)shaded exclusion region. 

D. Cold Dark Matter produced by misalignment 

As for axions produced by string decay, the relic den- 
sity of axions produced by misalignment should not ex- 
ceed the total cold dark matter density. From Eq. (fT4|) 
we have 

„..,3(e»,„)(i^)"° = 

2.8 X 10^(e|cD) ( ^ ) ^ ^cd,„, (35) 

where we neglected the anharmonic correction factor. 
This inequality provides a stringent upper limit on fa if 
(®qcd) order one. In particular, in the case of com- 
plete quantum diffusion during inflation, we have seen 
at the end of Sec. Ill Gl that (Oq^p) can be replaced by 
1.2 7r2/3 and 

/,<2.5xlO"Gev(^^)'^'. (36) 

This constraint is shown in Figs. [21 [3l as a dotted line 
(assuming tJcdm = 0.12), and excludes the light green 
region. 

In the absence of efficient quantum diffusion, 6qcd 
could take any nearly homogeneous value in our Universe: 



so, it is possible in principle to assume that (Bqqq) is 
extremely small (this coincidence can be justified by an- 
thropic considerations^), and to relax the bound on fa- 
However, the mean square cannot be fine-tuned to be 
smaller than the amplitude of quantum fluctuations at 
the end of inflation. Using Eq. (|15p. we see that (6qcd) 
can only vary within the range 



30 



/ Hi, 



< (0qcd) < 



TT 



(38) 



This gives a model-independent constraint 

'Wcdm 30 V fa 



Mnf < 2.5 X lO'^'GcV 



0.12 AiV 



10i2GeV 



which holds only in the region where 

2 



30 



inf 



< 



TT 

y 



(39) 



(40) 



otherwise it should be replaced by ([5^ . This bound 
excludes the dark green region in Figs. [21 [3] (assuming 
ojcdm = 0.12 from WMAP 3rd year results). 



E. Isocurvature modes 

As mentioned before, the axionic field induces an 
isocurvature component in the CMB anisotropics that 
must be considered when constraining the model. In this 
work, we assume that axions are the only source of isocur- 
vature modes. Taking expression ()25p for the isocurva- 
ture fraction a, replacing i?a by Wa/wcdm and using Eq. 
we obtain 



0.9 X lO^Cfe 



cdm 



Mp 

fa 



5/6 



(41) 



^ There has been plausible speculations that our presence in the 
universe may not be uncorrelated with the values of the funda- 
mental parameters in our theories. Such anthropic arguments 
normally arise in terms of conditional probability distributions 
of particular observables. In particular, the axion abundance 
is a natural parameter that could be bounded by those argu- 
ments, see e.g. Refs. [l^ . ISIl where it is suggested that the 
initial misalignment angle should be such that the main CDM 
component be axionic. In this case, one has a concrete prediction 
for i?a = ^a/^cdm = li 3.nd therefore the initial misalignment 
angle is directly related to the axion mass, see Eq. II35I I. 



nia = 9ficV (Oqcd) 



6/7 



(37) 



Having full diffusion, (©qcd) ~ 1.27r^/3, implies rria ~ 30 ^eV, 
just within reach of present axion dark matter experiments. On 
the other hand, we might happen to live in an unusual region 
of the universe with an extremely low value of (©qcq)) a large 
value of fa and still Ra = 1. 



Parameter 


Prior probability range 


Mean (±2cr) 




(0.016,0.030) 


0.023 ± 0.001 


I^cdm 


(0.08,0.16) 


0.123 ±0.008 


6 


(1.0,1.1) 


1.04 ±0.01 


T 


(0.01,0.2) 


0.09 ± 0.04 


"■ad 


(0.85,1.1) 


0.97 ±0.03 


As 


(2.7,4.5) 


3.1 ±0.1 


\a\ 


(-1,1) 


< 0.08 (95% c.l.) 



TABLE I: Prior probability ranges and means of the sampled 
parameters. 



We will now fit this model to the cosmological data in 
order to provide updated bounds on a and Wcdm- The 
analysis will also give constraints on the curvature power 
spectrum (\TZ(k)\'^), from which one can derive a relation 
between and Afjnf, using equation ([20]) . Therefore, 
our bounds on a and Wcdm will provide constraints in the 
(Minf, fa) plane. 

We assume a flat ACDM universe, with 3 species of 
massless neutrinos and seven free parameters, with the 
same notation as in Ref. [39| : uib, the physical baryon 
density, Wcdm, the total cold dark matter density (an ar- 
bitrary fraction of which is made of axions), 9, the angu- 
lar diameter of the sound horizon at decoupling, t, the 
optical depth to reionization, riad, the adiabatic spectral 
tilt. As = ln[10^"fc3(|7^p)], the overaU normalization of 
super-Hubble curvature perturbations during radiation 
domination at the pivot scale k = 0.002 Mpc^^ and a, 
the isocurvature contribution (we sampled from \a\ in- 
stead to avoid boundary effects near the maximum like- 
lihood region). Our data consists of CMB data (from 
WMAP (TT, TE and EE) [H, VSA [H, CBI ^ and 
ACBAR HBl); large scale structure data (2dFGRS JiS 
and the SDSS [1^); and supernovae data from Ref. [49j |. 
Note that previous studies indicated a very weak 
sensitivity of this data to niso, while in the present case 
?iiso = 1 — 2efc is very close to one, since the data favour 
Cfe <C 1. Thus, we safely fix niso to exactly one without 
modifying the results. Note also that the isocurvature 
mode arises exclusively from perturbations in the axion 
field while the adiabatic mode emerges from quantum 
fluctuations of the inflaton. Therefore, both contribu- 
tions are completely uncorrelated and we do not need to 
include an extra parameter that would measure this con- 
tribution (such as P under some parametrizations [38| ) 
in our analysis. 

A top-hat prior probability distribution was assigned 
to each parameter inside the ranges described in table [D 

We used the Metropolis-Hastings algorithm imple- 
mented by the publicly available code CosmoMC [501 to 
obtain 32 Monte Carlo Markov chains, getting a total of 
1.1 X 10^ samples. We find a x^/d.o.f. = 1.01 and the 
worst variance of chain means over the mean of chain 
variances value is 1.04 [sij . 

The one dimensional posterior probability distribu- 




tions for sampled and derived parameters are depicted 
in Fig. [T] In particular, the bestfit value for a is 
a = 8 X 10~^, and the 2a bound for the marginalized 
distribution is: 

a < 0.08 at 95% c.l (42) 

while for Wcdm we get: 

Wcdm = 0.123 ± 0.008 at 95% c.l. (43) 

Finally, our results for the amplitude of the primordial 
curvature spectrum gives the relation 



efc - 3 X 10* 



8 / Mnf 



Substituting in Eq. (|^T|) . we see that the bound a < 0.08 
finally provides the constraint 



Minf < lO^^GeV 



fa 



1012 GeV 



5/24 



corresponding to the central (yellow-) shaded forbidden 
region in Figs. O [31 As explained in Section 1111) this 
bound does not apply when thermal fluctuations induce 
PQ symmetry restoration after reheating, and we also 
assume -as in the rest of the literature- that it does 
not hold when quantum de Sitter fluctuations induce PQ 
symmetry restoration during inflation (although we be- 
lieve that this issue is not completely clear and deserves 
further study). So, for sufficiently large -ffjnf and/or ecs, 
the isocurvature constraint does not apply above a given 
line in the (Afjnf, fa) plane, and the allowed region is 
split in two parts (as in Fig. [3]). 



DISCUSSION 



General case 



The constraints discussed in the last section -and in 
particular the isocurvature mode limit- exclude a large 
region in parameter space, and preserve only two regions. 
The first one depends on the reheating efficiency'^, 

10^" GcV < fa < 1.2 X 10" GeV , 
5 X 10i'*F(eoff, fa) GeV < Mnf < 3 x 10^^ GeV , (46) 

with 



fa 



(44) V l0i"GeV 



'1.4x10-8 fa 



Eoff 1010 GeV 

The second one, for which the isocurvature mode is too 
small to be excluded by current cosmological data, cor- 
responds to 



(45) 10^" GeV < /a < 2.5x10" GeV, Mnf < 8x10^^ Qj,v , 

(48) 



where we assumed that the average misalignment angle 
in the observable universe is of order one: otherwise, the 
upper bound on fa could be relaxed significantly, while 
that on Minf would only increase slightly, as fa to the 
power 5/24, sec Eq. 

For a large class of well-known inflationary models, 
Afinf is typically of the order of 10^^ or lO^^ GeV, and the 



^ Here, we are still assuming for simplicity that A ~ 0.1. If this 
is not the case, the correct result is obtained by replacing the 
factor Cofi by eoff/{12A). 



relevant allowed window (if any) is the first one. In the 
next section, we feature the example of chaotic inflation 
with a quadratic potential. However, the second window 
PS)) is also relevant, since it is possible to build low- 
scale inflationary model compatible with WMAP con- 
straints and satisfying Afinf < 8x 10^^ GeV. Actually, the 
main motivation for the original hybrid inflation model 
of Ref. [ll] was precisely the possibility to evade axionic 
isocurvature constraints. Nowadays, many low scale in- 
flation models can be built in the generic framework of 
hybrid inflation (see Ref. [H^l for a review, or [s^ for a 
recent example). In this ease, the isocurvature perturba- 
tions carried by the axion are too small to be constrained 
by current data, and under the conditions of (jlSj) the PQ 
axion model is perfectly viable. 




fa (QeV) 



B. Bounds for chaotic inflation with a quadratic 
potential 

Inflation with a monomial potential V^(0) oc 0° is usu- 
ally called chaotic inflation. The latest WMAP results 
combined with other data sets essentially rule out cases 
with a > 4, while the quartic case a = 4 is only in 
marginal agreement with the data. Therefore, in this 
section we only consider the case of a quadratic potential 
V{(j)) = \vn?(jP, still favored by observational bounds. 

Using the COBE normalization, it is possible to prove 
that the mass should be of order m ^ 1.3 x 10~^ Mp, and 
to show that e-folds before the end of inflation, 

4>l = . (49) 

In particular, between A^obs and A^qcd, the scale of in- 
flation should be in the range 

1.5 X 10^6 GeV < Mnf < 2 x lO^*^ GeV , 
8 X 10^3 GeV < iJinf < 1 X 10^'' GeV . (50) 

So, if fa is close to or 10" GeV, one has i7i„f/27r > 
fa and the PQ symmetry is broken by quantum fluctu- 
ations during inflation. If these fluctuations erase the 
isocurvature modes on cosmological scales as suggested 
by Ref. [|, then chaotic inflation is compatible with the 
axion scenario. 

We could hope to find a second allowed window in the 
case where fa is very large and the misalignment angle 
is fine-tuned to a very small value (i.e., inside the light 
green region in Figs. [21 Note that this is not pos- 
sible, at least in the case of chaotic inflation. Indeed, 
reheating after chaotic inflation is expected to lead to a 
temperature Tj-h of the order of 10^° or 10^^ GeV. So, 
if fa is much bigger than 10^^ GeV, thermal symmetry 
restoration after inflation cannot take place, and isocur- 
vature bounds are applicable. Given the value of Afinf in 
Eq. ([50]). the isocurvature constraint ([45]) would enforce 
/a > 7 X GeV > M-p, which is not reahstic. 



FIG. 2: Bounds in the (Afinf, fa) plane, assuming that re- 
heating is very inefficient and the PQ symmetry can never be 
thermally restored (ecft < 4 x 10~^^). The energy scale of in- 
flation Afinf is bounded from above by r < 0.3 (top hatched) 
and from below by the requirement of successful baryogenesis 
(bottom hatched). Possible values of the PQ scale /„ lie be- 
tween the regions excluded by supernovae (left blue) and by 
the bound on the axionic relic density, produced by one of the 
following two mechanisms: axionic string decay (red) or mis- 
alignment angle at the time of the QCD transition (green). 
For the bounds related to the relic density, we use a solid line 
for the model independent bound (with a minimal, fine-tuned 
value of the misalignment angle), and a dotted line for the 
model dependent bound (with a generic misalignment angle 
of order one). The yellow region is excluded by the limit on 
the isocurvature mode amplitude a: this is the main result of 
this work. If large quantum fluctuations during inflation can 
erase isocurvature modes on cosmological scale (as suggested 
by Ref. 0]), the isocurvature bound does not apply above the 
line corresponding to fa = ffinf/Svr. The remaining allowed 
regions are left in white. Details are explained in the text. 



C. Possible loopholes 

In this subsection we will explore those loopholes we 
have left open for the axion to be the dominant compo- 
nent of cold dark matter. 



1. Production of cosmic strings during reheating 

Even if the reheating temperature of the universe is too 
low for a thermal phase transition at the Peccei-Quinn 
scale, it is possible that axionic cosmic strings be formed 
during preheating if the field responsible for symmetry 
breaking at the end of inflation is the Peceei-Quinn field. 
Then, the residual global U{1) symmetry of the vacuum 
gives rise to axionic cosmic strings. At present there 
is no prediction for what is the scaling limit of such a 
mesh of strings produced during preheating. A crucial 
quantity that requires evaluation is the fraction of energy 



m_, (eV) 

10"^ 10"* 10"^ 




f, (QeV) 



FIG. 3: Same as previous figure, but assuming now that 
reheating is very efficient (teff ~ 10~*) so tliat the the PQ 
symmetry is thermally restored above the dotted line. In 
this case, the constraint from string production and decay 
(red) extends to smaller values of Mini, but above the dotted 
line and for 10^°GeV < fa < 1.2 x 10"GeV there is a new 
allowed region (corresponding to the erasure of isocurvature 
perturbations in the thermal bath). 



density in infinite strings produced at preheating, since 
they are the ones that will give the largest contribution to 
the axion energy density. It could then be that axionic 
strings produced at preheating may be responsible for 
the present axion abundance. We leave for the future 
the investigation of this interesting possibility. 



2. Axion dilution by late inflation 

When imposing bounds on the axion mass from its 
present abundance, it is assumed that no there is no sig- 
nificant late entropy production [s^ ] or dilution from a 
secondary stage of inflation. We truly don't know. It 
has been speculated that a short period of inflation may 
be required for clcctroweak baryogenesis to proceed [55| . 
In such a case, a few c-folds (N ~ 5) of late inflation 
may dilute the actual axion energy density by a factor 
ry _ g-3N ^ IQ-'', easily evading the bounds. 



VI. CONCLUSIONS 

In this paper, we reanalyzed the bounds on the cos- 
mological scenario in which the cold dark matter is com- 
posed of axions plus some other component (like e.g. ncu- 
tralinos), and the energy density of axions is produced 
by the misalignment mechanism at the time of the QCD 
transition. In this fraction of the cosmological 

perturbations consists of isocurvature modes related to 
the quantum perturbations of the axion during inflation. 
We use this possible signature plus the contribution to 
dark matter energy density as a tracer of axions in the 
universe. In that sense, this work is similar to that pre- 
sented in (56j . However, we make a stronger statement 
about the generation of isocurvature modes in the case in 
which the PQ symmetry is restored by quantum fluctua- 
tions during inflation, and we signiflcantly improve pre- 
vious bounds coming from the non-observation of isocur- 
vature modes in the CMB, in particular, given the recent 
WMAP 3-year data. 

A narrow window for the PQ scale 10^° GeV < fa < 
10^^ GeV still remains open, but we show that the main 
consequence of recent data on cosmological perturba- 
tions is to limit the possible energy scale of inflation in 
this context, in order not to have an excessively large 
contribution of isocurvature perturbations to the CMB 
anisotropics. In particular, we show that in the axion 
scenario, the energy scale of inflation cannot be in the 
range 8 x 10^^ < A'/inf < 5 x lO" GeV, unless one 
of the two situations occurs: either reheating leads to 
a temperature T^h > fa, with an efficiency parameter 
Ecff > 4 x 10~^^, and there is another allowed region 
with 6 X lO^^e^s^^ GeV < Afi^f < 3 x 10^^ GeV; or 
the misalignment angle is fine-tuned to very small values 
(this coincidence can be motivated by anthropic consid- 
erations) and the upper bound Minf < 10^"^ GeV can be 
slightly weakened (by at most one order of magnitude) . 

These bounds may be of interest taking into account 
particle physics experiments searching for the axion, 
which may help to put very stringent constraints on in- 
flationary models. Detecting the QCD axion coifld shed 
some light on the scale of inflation and possibly into the 
mechanism responsible for inflation. 

In this respect, there is an intriguing possibility that 
the PVLAS experiment [13] may have observed a pseudo- 
scalar particle coupled to photons. The nature of this 
particle is yet to be decided, since its properties seem 
in conflict with present bounds on the axion coupling to 
matter [l^, [13] , see however [H, [s^ . 



3. Coupling between radial part of PQ field and other fields. 

The radial part of the PQ field could interact with 
other fields, including the inflaton (see e.g. [l^). Then 
the effective fa during inflation might be anything, and 
the isocurvature bounds might be alleviated. 
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